### Subset for survey
dsub <- subset(srs_design, M_20 == 1)

### 
svyby(~age, ~M_20, srs_design, svymean)

### P value - cuts 
svyttest(I5_P11_1 ~ wariant, dsub)

svyttest(I5_P11_2 ~ wariant, dsub)

svyttest(I5_P11_3 ~ wariant, dsub)

svyttest(I5_P11_4 ~ wariant, dsub)

svyttest(I5_P11_5 ~ wariant, dsub)

svyttest(I5_P11_6 ~ wariant, dsub)

svyttest(I5_P11_7 ~ wariant, dsub)

svyttest(I5_P11_8 ~ wariant, dsub)

### P value - hikes 

svyttest(I6_P12_1 ~ wariant, dsub)

svyttest(I6_P12_2 ~ wariant, dsub)

svyttest(I6_P12_3 ~ wariant, dsub)

 svyttest(I6_P12_4 ~ wariant, dsub)

svyttest(I6_P12_5 ~ wariant, dsub)

svyttest(I6_P12_6 ~ wariant, dsub)

svyttest(I6_P12_7 ~ wariant, dsub)

svyttest(I6_P12_8 ~ wariant, dsub)
